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ABSTRACT 

Directional melting of binary systems, as encountered during 
seeding in melt growth, is analyzed for concurrent compositional 
changes at the crystal-melt interface. It is shown that steady 
state conditions cannot normally be reached during seeding and 
that the growth interface temperature at the initial stages of 
seeded growth is a function of backmelt conditions. The theoretical 
treatment is numerically applied to Hg. Cd Te and Ga-doped Ge. 
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1. INTRODUCTION 


Directional melting, contrary to solidification, has been the subject of 
only very few studies. Woodruff and Forty (1) investigated melting and 
freezing of a Na-K alloy by in situ transmission UV microscopy; optical 
phenomena observed in their melting experiments have been interpreted as 

I 

indicative for Na diffusion from the solid into the adjacent liquid phase. The 
directional melting studies by Verhoeven et al. (2) in Sn-Bi, carried out at a 
rate of 0.2 pm/sec, reportedly revealed & composition gradient region of 650 pm 
in the solid ahead of the receding interface. The establishment of such an 
extended gradient region (boundary layer) was found by them to conflict with 
theoretical prediction (1,3 pm) based on the application of their steady state 
segregation equation (3) for melting. It should be noted that the extended 
solute depleted region in their experiments was a two-phase region with liquid 
droplets appearing dispersed in the solid phase. They also report that no 
solute depletion region in the solid could be observed at higher melting rates 
(20 pm/sec) . Verhoeven and Heimes (4) provided analytical solutions for 
initial melting transients with preceding solidification which consider 
diffusion in both the solid and liquid phases. 

All seeded crystal growth from the melt involves as a first step a varying 
degree of backmelting of a crystal. This approach, generally successful in 
systems with small liquidus-solidus separation such as doped semiconductors, 
must, however, be expected to be problematic if binary or multicomponent systems 
with large liquidus-solidus separation (Hg, Cd Te, Pb. Sn 'Te, Ge-Si and the 
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like) are involved. Any changes in composition that occur during melting at 
the crystal-melt interface are of importance since they control the growth 
interface temperature and its stability. Moreover, both parameters influence 
the system’s behavior during equilibration (melting arrest) and may affect the 


t 


degree of crystalline and compositional perfection achievable during ensuing 
solidification. Seeded crystal growth, for exmple by the Bridgman and 
Czochralski techniques, involves directional melting under boundary and 
initial conditions which were not considered in the earlier treatments (4). 

It is the purpose of this work to analyze the evolution of the solute 

I 

concentrations, c!(tu and C*(t), at a crystal-melt interface and the related 
interface temperature during melting, with emphasis on the initial transient 
behavior. The analysis applies to binary systems in a configuration where a 
crystal of composition Cg in contact with a melt of composition C L is direc- 

i 

tionally melted. The theoretical treatment is carried out in non-directional 
form; it is numerically applied to Hg Jx Cd x Te and Ga-doped Ge to provide for 
comparative analysis of systems with small and large liquidus-solidus separa- 
tion. 

2. Theoretical Analysis of Solute Redistribution During Directional Melting 

2 . 1 Definition of boundary conditions and mass transport relationships 

The analysis for directional .melting of binary solids is made with the 
following assumptions: 

(a) the crystal-melt interface is planar; one-dimensional 
analysis; 

♦ ■ 

(b) the interface displacement R is constant; 

(O the mass transports in the liquid and in the solid 
obey Pick's law; 

(d) convection is neglected; 

(e) the interface distribution coefficient k, taken to 
be identical with the equilibrium distribution 
coefficient, is a constant; 

(f) the diffusion coefficients in the liquid (D^) and 
in the solid (Dg) are constant; 
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(g) the initial length of the solid is larger than the 
characteristic length of diffusion; 

(h) the initial compositions of the solid and contacting 
melt are in thermodynamic equilibrium (C s * kC L ). 

• % 

It is recognized that a constant D is not realized during melting since 

steep temperature and concentration gradients encountered across the crystal - 
molt interface must affect its magnitude. A detailed determination of the 
changing driving forces resulting from the existing gradients is complex 
because of their interdependence. Presently, the contributions of these 
driving forces are accounted for as an additional flux term, Jp(z,t), which 
modifies Pick’s second law as follows: 

V t • C(z,t) * -V z • J(z,t) 

with 

J(z,t) = -DV • C(z,t) - <D b >V, • C(z,t) 

Z r Z 

or 

J(z,t) = -D eff V z • C(z,t) 

where D e ££ « d + <D p > is defined as an "effective" coefficient of diffusion. 
The concept of D e ££ has been used earlier (5) to account for the effect of 
convection. In this paper the notation Dg is used for D e ££ in the solid and 
for D e ££ in the melt, With regard to assumption (e) , it is customary in 
dilute binary systems such as doped semiconductors with small liquidus-solidus 
separation to take k to be independent of concentration. This approximation 
is in principle erroneous for systems with large liquidus-solidus separation. 
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The concentration dependence of k is known from the phase diagram and it could 
be added to the set of boundary conditions (providing that the displacement 
rate of the interface is constant • see Appendix). Since no general analytical 
treatment of the relationship between k(c,t) and R(t) is as yet available, R 

and k are taken to be constant, 

» % 

2.2 Development of composition transients associated with directional melting 
Making use of the preceding assumptions, Fick’s second law, formulated for 
a coordinate system which moves with the interface according to the configura- 
tion shown in Fig. 1, assumes the form: 


3 2 C s dC s 

3C s 
' at 

[1] 

3 2 C. 3C. 

3C, 



_ L 

* *3T 

t 

[2] 

and boundary conditions 

for equations [1] and [2] are: 


t * 0 C c = C 

a 0 

ze]0,+«] 

[3] 

C. * C /k 
L 0 

Z€]0,-«>] 

[4] 

interface (z * 0) : 



C S ■ k C *L 

¥• t > 0 

[5] 


Solute flux conservation through the interface yields: 


3C, 


J S 3z 


3C. 


z=0 + 


♦ R c s ■ d l TT 


+ r c; v - 1 > o 


z=0" 


[ 6 ] 


Using the change of scale: 


C «= C - C 


[7] 
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and the dimensionless variables: 


K w 


T ■ 4-t [9] 

U L 


« - 5^ UO] 
U L 


the system of eqs. [1] and [2] becomes: 


& C S S 3c s k 
° 3y 2 + ly“ * 9 t 


[ 11 ] 


i\ 3C l 3C l 


9y 


r + -sr* it 


[12] 


with the corresponding initial and boundary conditions: 


T ■ 0 

C s • 0 

yejo. 400 ] 

[13] 
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ye]0,-«] 

[14] 

the interface (y * 0) : 



Cj * kCf + C (k - 1) 

w L 0 

V- T > 0 

[15] 

SC c 

| 



s 

ay 

+ - Li 

y=0* 5 * |y= 

+ cf . 
0‘ L 

[16] 


The set of equations [11-12] with the boundary conditions [14 to 16] and 
the initial conditions [13] can be solved by using the Laplace transforms: 


. W’ 


C s e -*dx 


[17] 


and 


Yr = 


r 


C L e” ST dT 


[18] 
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yielding: 



fs 

dy 



0 


d\ d *L 


with 


q - C„(i - 1) 

Taking the general solutions: 

Y S ■ A o-^ /a 


119) 


[ 20 ] 


[ 22 ) 

[23] 


and using [13 to 16], one obtains: 



kq expC-l-^X 


[24] 


and 



+ 


k 


exp(-~ + 


JL 1 k 

2a ” 2 ” "572 
a 


t+*)y 

y F*> -i/F* 


) 


[25] 


These relationships depend mainly on the two parameters k and a = Dg/D^. The 
inverse Laplace transforms are obtained in the general case by a numerical 
method. Analytical solutions will be given for particular values of a. 

The inverse Laplace transform is made by approximating the integral: 
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C(y,T) 


l 

3iT 


/c+i» 

c-i« 


e 4T ?(4)d4 


[26) 


by a trapezoidal rule using discrete values of the Laplace transform ?(4) along 
the integration path A * c ♦ iw. For this purpose the integral [26] is approx- 
imated by Fourier cosine series which are evaluated using the Fast Fourier 

. . % 

transform algorithm. The numerical inversion of [24] gives C^(y,x) for 
given values of o and k. Figure 2 is a plot of t q gg , the time necessary to 
reach 99* of steady state as a function of k, for four different values of a 
obtained by the inversion of [24]. [The real time (t) is related to the dimen- 
sionless time x through t * Dji/R ,] 

At the hypothetical upper limit (a c 1), the liquid and the solid behave 
symmetrically and a corresponding solution of [24] is given by: 


c s (y,T) » ~(1 - k)erfc(^p -~)exp(-y) - erfc(^p- ~) [27] 


/F, 


which reduces at the interface to: 


C£(0,x) C Q (k - 1) erf ~ 

Equation [28 j confirms the plotted data, 99% of steady state (erf 


/F 


[28] 


0.99) is 


reached for all values of Ml, with t = 13.4. The corresponding curve on 
Fig. 2 is used as a reference line to determine the variation of x g gg with a 
for any given value of k. 

Considering the special case of a -► 0, the inversion of [24] yields: 


c s (y,T) = {(1 -k)erfc(^-^)exp(~y) + £erfc 


/F, 


- T)erfc( 2^F + iir 1293 

k 


#• 
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PAqf 

cp P0 °« «S£ 


C^(0,t) ■ y (Cl *k)wfc(-^S ^ ^crfc (y-) 


- (^expC^* T)erfc(^|^ t) 


[30J 


2.3 Discussion of segregation analysis 

The relationships [29] and [30] can be applied to Solidification by using 

c c c 

the transform S defined by k — ► 1/k, S — * L, L — ♦> S. It yields a solution 
which is identical with that given by Smith et al. (7) where D g is considered 
negligibly small. Figure 2 shows that! curves drawn for different values of 

a * 

a convorge to equation [29] for a 5 10 . Generally accepted values for D g and 

i a R 7 A 

in semiconductors lie in the range 10 - 10 cm /sec for D g and 5x10 - 
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10 cm /sec for D. . l^ci/cfore the Smith et al. solution (or its equivalent for 

h 

melting [29]) must be considered a good approximation. It should be pointed 
out, however, that as yet there are only limited data (of reasonable precision) 
for diffusion coefficients at melting point temperatures. While the uncertainty 
concerning the magnitude of is not expected to have a noticeable effect on 
the value of a, D g can be considerably larger than determined by extrapolation 
from , lower temperatures. The temperature controlled increase in vacancy con- 
centration, for example, will affect D c and may thus have a significant effect 

M 

on the value of a. The data show that already for a * 10 the Smith et al. 
approximation (or eq. [29]) no longer applies, and that for any given interface 
displacement rate, the compositional transients in the solid (Tq g g) (and in 
the liquid) increase with increasing a for k < 2 and decrease for k > 2. 

2.4 Implications of segregation behavior on seeded crystal growth 

* 

It is of interest to consider the present analysis in terms of real time 
(t) behavior, i.c. to consider, for example, the specific effect of the rate of 
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directional molting (R) on the time required to reach steady state. For this 
purpose, it ig convenient to analyze the behavior of specific systems. An 
arbitrary choice is Hg^Cd^Te of the composition x * 0.2 (C Q * 0,2). Figure 
3 presents th« concentration of CdTe in the solid (HgCdTe) at the crystal-melt 
interface (Cj) as a function of t for k g and k^. In real time (t), the CdTe 
concentration (C g ) is dependent on the rate of melting (R) which, in order to 
be meaningful, must satisfy the stability criteria (3). At the beginning of 
melting, the system is at equilibrium [C c (z,0) * C , C.(z,0) * C /k] , and 

therefore under conditions of absolute stability. As melting proceeds, the 

» 

interface stability is affected by out-diffusion and the subsequent formation 
of the boundary layer in the solid. Thus, while for C, (Q,z) * C /k melting 
can be initiated at any rate, out -diffusion and interface breakdown must be 
anticipated if the rate exceeds a critical value. If the stability criterion 
is strictly applied to a steady state melting configuration, the G g /R ratio 

g 

required for stability is very high («10 °C/sec for Hg () g C4 0 2 Te). Suctl a high ' 
value is achieved by either a very high thermal gradient or a very slow rate 
of interface displacement. Making use of Fig, 3 and applying numerical values 
for HgCdTe, the present analysis shows that under neither condition will 
steady state be reached during seeding in Bridgman or Czochralski 
growth whore the seed is back-melted for several millimeters only (see 

Table I); the crystal-melt interface temperature during molting must therefore be 
expected to lie somewhere between the solidus and the liquidus temperatures. 

Considering the back-molting of an ingot of Hg fl Cd n ^Te for half an hour at 

-2 

the rate of 3,7 x 10 pm/ sec (the rate required for absolute stability at 

steady state if G g * 1000°C/cm) , the x value at the end of back-melting (i.e. 
after half an hour) is 4 x 10" 4 and, according to Fig, 3, the concentration in 
the solid (C g ) will bo essentially unchanged at C Q and the interface temperature 
will remain at about T g . 


Hg t - x Cd 

* k, = 2.6 


(3 0 ) 9in|DJ0dujax 
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The dependence of the length of transient behavior on the magnitude of 
the liquidus-solidus separation can be demonstrated by comparing Hg Q gCd ( . Je 
with Ga-doped Ge where the liquidus-solidus separation is very 

small (less than 1 # C). The x curves of Ga-doped Ge are presented in Fig. 4. 

18 3 

At the composition 10 at/cm of Ga, the G g /R ratio assumes the value of 

2.9 x 10^ # C/sec (m g • 0.366 # C/%; D g * 10"® cm^/sec; k » 0.087), Since the 

Ga-doped Ge system has been extensively studied, the credibility of the data 

used for the different parameters is much higher than for HgCdTe. Table I 

shows that the length to be melted to reach steady state (z m ) in doped Ge is 

much smaller than in Hg_ e Cd rt «Te. The interface concentration in the solid 

under the conditions depicted in Table I can be calculated using the x curves 

in Fig. 4. Assuming again that back-melting is to proceed for 0.5 hours at 

G g «* 10°C/cm and R = 34 pm/sec, we obtain a x value of 99 which indicates 

that steady state has been reached and that the concentration in the solid 

at the interface, C g ,,is kC Q . For this system, x Q (steady state) is already 

reached after about 9 minutes. With continuing melting, the crystal-melt 

interface temperature changed from T g to T^, the liquidus temperature. It is 

noteworthy that for systems with a liquidus-solidus separation smaller than 

in Hg , Q Cd n 0 Te but larger than in very dilute systems, the temperature of the 

interface (T*) at the end of back-melting is expected to be somewhere 

between T g and T^J its value, at a given time, can be calculated from 

appropriate t curves. The interface temperature is related to the 

composition at the interface through the relation T* s T m + ro t G g , where T^ is 

the melting temperature of the major constituent (for example, HgTe in 

Hg n 0 Cd n JTe) and m* is the slope of the straight line joining T and the point 
u t o u«2 i in 

on the solidus curve corresponding to the concentration C g on the phase 
diagram. Therefore, the deviation of the interface temperature from the 
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solidus temperature during transient melting (AT* * T* - T s ) in non-dimensional 
notation is given by: 


AT* « m*c£ - m, 


t^S '"t-0 o 


[31] 


As Cg approaches the value of kC Ql AT* approaches a maximum value for 
which T* * T l . ‘ % ’ 

The preceding analysis and conclusions apply to initial conditions whereby 
a (binary) solid is in thermodynamic equilibrium with the contacting melt 
(Cg * kC L ). This condition can in principle be realized during seeding in 
Czochralski and Bridgman-type growth. The analysis as such, however, is not 
applicable to configurations in which initially the entire charge is solid and 
seeded "re-growth" is accomplished by first directionally melting a fraction 
of the solid [configuration used during growth experiments in reduced gravity 
environment (8)) . Initial molting in this case cannot a priori be expected to 
take place at T g . If the phase transformation is initiated under diffusionless 
conditions, melting will take place at T g < Tq < T^, a temperature at which the 
solid and the liquid formed have equal free energies (9) but are at a different 
chemical potential. Maintaining this kinotically controlled metastabl 3 state 
will obviously bo a function of the experimental conditions under which melting 
proceeds. Depending on the extent of diffusional adjustment of the phase 
boundary composition toward equilibrium, a transient in the crystal -melt inter- 
face temperature will be encountered. These transients to equilibrium for 
systems with t C g /k arc currently under investigation. 


APPENDIX 

As a consequence of segregation during transient molting or solidification 
the growth-interface temperature (T*) is subject to change and the interface 


